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We present some recent developments in our study of inhomogeneous 
chiral symmetry breaking phases in the Nambu- Jona-Lasinio model. First, 
we investigate different kinds of one- and two-dimensional spatial modu- 
lations of the chiral condensate within the inhomogeneous "island" and 
compare their free energies. Next, we employ the Polyakov-loop extended 
version of the model to study the effects of varying the number of colors 
on the inhomogeneous region. Finally, we discuss the properties of an in- 
homogeneous "continent" which appears in our model at higher chemical 
potentials, and analyze its origin. 

PACS numbers: 21.65. Qr,12.38.Mh,12.39.Fe 

1. Introduction 

The study of spatially modulated ground states in strongly interacting 
systems is an old topic [H [21 [3l H] (see Ref. for a recent review) which has 
recently received new attention. It has been argued some time ago that the 
favored ground state of a dense Fermi sea of quarks should be characterized 
by a spatial modulation of the chiral condensate, at least in the limit of a 
large number of colors (Nc) [SI [Z|- More recent studies, especially in the 
context of quarkyonic matter, seem to support this hypothesis [HI [9| [TO]. 

But also for physical case of three colors, Nambu- Jona-Lasinio- (NJL-) 
type model studies have revealed the presence of an inhomogeneous "island" 
at intermediate chemical potentials and low temperatures, namely in the 
region where the usual first-order chiral phase transition would occur when 
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limiting to homogeneous phases only |1H I12 [ [T3]. In particular, the critical 
endpoint of the phase boundary is covered by the inhomogeneous region 
and thus disappears from the phase diagram |14j . 

In these proceedings we present some of the most recent results in our 
NJL- model study of inhomogeneous chiral symmetry breaking phases. After 
briefly outlining our methods for tackling this problem and the technical 
difficulties associated with it, we show numerical results for different kinds 
of one- and two-dimensional modulations and compare their free energies. 
In order to build a bridge towards large- A'^c studies, we then consider a 
Polyakov-loop extended NJL model for an arbitrary (large) number of colors 
and see how this modification alters the inhomogeneous phase. Finally, we 
discuss the properties of a second inhomogeneous region, which appears in 
our model at higher chemical potentials, and try to understand the origin 
of this new "continent". In particular, we investigate whether it is to be 
interpreted as a regularization artifact or as a simple consequence of the 
model interaction. 



Starting from the two-flavor Nambu-Jona Lasinio Lagrangian |15j . 



we perform the mean-field approximation, allowing for a spatial dependence 
of the scalar and pseudoscalar condensates: 



The mean-field Lagrangian can be rewritten by introducing an effective 
Hamilton operator % [131 US] • 



2. Inhomogeneous phases in NJL 



(1) 



(^V> = S{x) 



(2) 



(3) 



with T-L defined as 



n = -f^ ij-d + m- 2G{S + i-f^TaP) 



(4) 



Explicitly, we can write down T-L as [13) 





where the chiral representation for the Dirac matrices has been used and 
we introduced an effective inhomogeneous "mass" function 
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M{x) = m- 2G{S{x) + iP{x)) . (6) 

In order to evaluate the thermodynamic potential associated to this (so far 
generic) spatially modulated chiral condensate, we work, as customary, in 
imaginary time and switch to momentum space. By assuming static (i.e. 
time-independent) condensates, we are able to perform explicitly the sum 
over Matsubara frequencies and obtain the expression (up to a constant) 

^t{T,^i■M{x)) = ^kva{T,IJi\M{x)) + rt,ond{M{x)), (7) 

with 

^cona{M{x)) = 1 1^ d^X "^1' , (8) 

where V is the volume of the system, and 

^kvn{T,ii;M{x)) = -T^log (^2 cosh (^^|^)) , (9) 

En 

where the sum is over all eigenvalues of % in color, flavor, Dirac and 
momentum space. 

For homogeneous conde nsates the diagonalization of % in Dirac space 
simply gives twice Ep = ± + and the eigenvalue sum may be triv- 
ially turned into an integral over all possible momenta. In presence of an 
inhomogcneous condensate, however, the diagonalization of H is a highly 
non-trivial task, since the quarks may exchange momenta by scattering off 
the inhomogcneous condensate. This means that the spatially modulated 
mass term will effectively couple quarks with different momenta and the 
resulting structure will not be diagonal in momentum space. 

In the following we assume to have a lattice structure and, thus, a peri- 
odic shape of the modulation. This means that we can expand the modu- 
lated chiral condensate in a Fourier series, 

M{x) = Y,M^,e^^>'-\ (10) 

with discrete momenta g^, which form a reciprocal lattice (RL). A generic 
element of T-L in momentum space then takes the form 

nj _ ( ~^ ' Pm^pm^Pn '}2iqf,^qk^Pn,Pm+qk \ /-i -i \ 

\ Z^qk ^Qk Pn,Pm-qk " yrn"pm.,Pn J 

where runs over all momenta of the RL, making obvious the non- 
diagonal structure of the matrix. In turn, momenta which do not differ by 
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an element of the RL are not coupled, so that % can be decomposed into 
a block diagonal form, where each block can be labelled by a momentum 
of the first Brillouin zone (BZ). This allows to decompose the eigenvalue 
sum in Eq. d?]) into a momentum integral over the BZ times a sum over the 
discrete eigenvalues of each block. 

For the simplest case of the so-called chiral density wave (CDW), namely 
a plane wave characterized by a single momenturr0 

M{x) = Me'^-^, (12) 

the sum in ()lip is then obviously given by a single term, and the resulting 
BZ-projected matrix within this framework is characterized by a narrow 
band structure, with only the first off-diagonal block filled with nonzero 
entries. 

For general periodic structures, although the numerical diagonalization 
procedure is in principle straightforward, its practical implementation turns 
out to be computationally demanding. Therefore, in order to simplify the 
problem, we limit the generality of our ansatz Eq. (jlOp for the spatially mod- 
ulated condensate and consider simpler, lower-dimensional modulations. 

3. One-dimensional modulations 

The eigenvalue problem simplifies considerably when the chiral con- 
densate is allowed to vary only in one spatial dimension. In this special 
case it has been observed that the dimensionally reduced effective Hamil- 
tonian becomes formally identical to that of the 1+1-dimensional Gross- 
Neveu model [13] , for which self-consistent solutions are already well known 

[T8| [T9l [20] . One finds (from now on we will consider only results in the 
chiral limit for simplicity) 

M{z) = V^Asn(Az|i/) , (13) 

where sn{^z\v) is a Jacobi elliptic function |21j . For this kind of solutions, 
an analytical expression for the eigenvalue spectrum, depending on the ellip- 
tic parameter v and the amplitude of the condensate A can be obtained and 
the minimization of the thermodynamic potential may easily be performed 
with respect to these two quantities [131. No numerical diagonalization of 
% is thus needed. In particular, Eq. ([9|) can be written as 

^ Some authors introduce an additional factor of two in the exponent, M{x) — Me*'^''^. 
This is motivated by the fact that the favored value of |(5| is roughly of the order 
of 2/Lj, so that \q\ is of the order of /x. Here we prefer the definition without the 
factor of two, which is more consistent with the general ansatz of Eq. (|1UI) and other 
modulations studied in this article. 
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Fig. 1. Left: Phase diagram obtained when allowing for solitonic solutions. The 
black dotted lines represent the second-order transition lines separating the ho- 
mogeneous chirally broken phase (left), the chirally restored phase (right) and the 
inhomogeneous phase (shaded region). The blue solid line represents the first-order 
phase boundary obtained when limiting to homogeneous order parameters. It is 
completely covered by the inhomogeneous phase and ends at the "Lifshitz point" 
where the three second-order phase boundaries meet. Right: Mass function M{z) 
at T = for two different values of the chemical potential. 



f'OO 

^kvn = -NfN, / dEp{E){f,ac + fraed) , (14) 

Jo 

where all functions which enter the integrand are known analytically [13] : 
p{E) is the spectral density associated to this kind of solutions, f^ac is 
the (divergent) vacuum contribution related to the Dirac sea, which we 
regularize by introducing Pauli-Villars-type counterterms [22], 

3 

f,ac(.E,A) = J2cjVE^+j^\ Co = -C3 = l, C2 = -Ci = 3, (15) 
j=0 

and fmed describes the medium part: 

f^UE) = Tlog (l + e-(^-^)/^) + Tlog (l + e-(^+'^)/^) . (16) 

The resulting phase diagram is shown in Fig. [TJ It features a region at 
low temperatures and intermediate chemical potentials where an inhomoge- 
neous phase is favored over the homogeneous chirally broken and restored 
solutions. At the onset of this inhomogeneous island the chiral condensate 
assumes the shape of a single soliton (which is thermodynamically degen- 
erate with the homogeneous broken solution) and then progressively varies 
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Fig. 2. Amplitude and wave number at T = as functions of the chemical potential 
for the CDW ansatz, Eq. (jl2l) (left) and for sinusoidal modulations, Eq. (IT7|) (right). 



into a more sinusoidal shape until it gradually melts when reaching the 
restored phase. In this case, all phase transitions are second order [13\ I14j. 

It is of course possible to restrict the one-dimensional ansatz to some- 
thing even simpler, such as a chiral density wave ()12p |1H 112] or a real 
sinusoidal modulation, 

M{z) = Mcos{Qz) . (17) 

When implementing the CDW modulation one finds that the BZ-projected 
H is immediately block-diagonal and an analytical expression for the density 
of states can be obtained straightforwardly. The thermodynamic potential 
can then again be worked out using Eq. ()14p . with the only difference being 
a different expression for p{E) [13j. For the real sinusoidal modulation, on 
the other hand, a brute- force numerical diagonalization in Dirac and mo- 
mentum space is required. Results for the order parameters of these two 
kinds of modulation are shown in Fig. [21 After comparing the free energy 
of these modulations (see also section 14. 2p , one finds that both are favored 
over the homogeneous solutions in a very similar window in chemical poten- 
tial. In fact, it has been argued [23] that all these inhomogeneous solutions 
should share a common second order transition to the restored phase, while 
the onset from the chirally broken phase need not be the same, since in 
this case the two phases are separated by a first order lin^. Nevertheless, 
numerical results show how the onsets of the two phases occur at roughly 
the same value of chemical potential, with a slightly larger window for the 
real sinusoidal modulation. 



^ While the solitonic ansatz may assume a shape that is thermodynamically degenerate 
with the homogeneous broken case, there is no smooth way for these less general 
solutions to approach a spatially constant shape unless for the trivial case M — 0. 
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4. 2D modulations 

While the study of one-dimensional modulations has been able to pro- 
vide a first insight on the importance of spatial modulations of the chiral 
condensate, it is natural to expect that also higher-dimensional modulations 
could appear in the phase diagram of a 3+1 dimensional system. Aside from 
being a more general ansatz, higher-dimensional modulations are also un- 
affected by the instability with respect to fluctuations which prevents the 
formation of a true ID crystalline structure at finite temperature [23] . 

The procedure for handling a higher-dimensional modulation is in princi- 
ple identical to that already outlined for the one-dimensional case. However, 
for higher dimensions no straightforward analytical results are available and 
the numerical diagonalization involves a bigger matrix, since one is dealing 
with two- or three-dimensional momenta in %. 

An obvious next step is to study two-dimensional modulations. Then 
the thermodynamic potential takes the form 



NfN, 



2 cosh 



2T 



+ n 



cond 



(18) 

where p± is the momentum perpendicular to the directions of the spatial 
modulation and one introduces E± = sgn{\£j ^ + p^, with being the 

eigenvalues of the dimensionally-reduced depending on the momentum 
G BZ. 



4-.1. Square modulation: egg carton 

Since there are no known results to guide the choice for the shape of the 
spatial modulation, we should in principle study different geometries of the 
two-dimensional crystal and for each of them minimize the thermodynamic 
potential with respect to the Fourier components on the corresponding 
RL. As a first step in this direction we consider a square lattice and a simple 
product of cosines. 



M(x, y) = A cos(Qx) cos{Qy) 



(19) 



thus restricting ourselves to the first harmonics only in each direction. This 
highly symmetric "egg carton" (shown in Fig. [3]) is possibly the simplest 
real modulation one could think of in two dimensions f| 



This ansatz is equivalent to a modulation of the form M{x',y') — A'[cos(Q'a;') + 
cos(Q'j/')], which is one of the shapes discussed in [T^. Here A' — A/2, Q' = V^Q 
and (x ,y') is related to {x,y) by a rotation of 7r/4 about the z-axis. 
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Fig. 3. Left: Shape of the egg-carton modulation. Right: Results for Q and A 
obtained by numerical minimization of the thermodynamic potential for this two- 
dimensional ansatz. 

After diagonalizing % numerically, the thermodynamic potential p8p 
is minimized with respect to the two parameters A and Q, characteriz- 
ing amplitude and period of the modulation, respectively. The results of 
the numerical minimization are presented on the right-hand side of Fig. [3l 
Their main features are in qualitative agreement with the one-dimensional 
examples discussed before. We find again a sharp onset around 310 
MeV and a smooth approach to the restored phase, which is reached at the 
same chemical potential as for the ID modulations via a second-order phase 
transition. 



The results presented in Figs. [2] and E] have been obtained by assuming 
a single ansatz for the mass modulation (CDW, cosine or egg carton) in 
each case. Under this restriction we found that the different inhomogeneous 
phases are energetically favored over the homogeneous solutions in a win- 
dow, which at T = and for our model parameters lies roughly between 
/X w 300 and 350 MeV. 

The next obvious step is now to compare the values of the thermody- 
namic potential of these solutions with each other in order to find out which 
of them has the lowest free energy. In fact, for the one-dimensional modula- 
tions, we know already that the CDW and the cosine are disfavored against 
the solitonic solutions, but there could still be higher-dimensional modula- 
tions, which have an even lower free energy than the latter. In particular, 
it has been argued in the context of quarkyonic matter studies that with 
increasing density such higher-dimensional solutions will be favored, at least 
in the limit of a large number of colors [3 [TU]. On the other hand, using 
Ginzburg-Landau (GL) arguments, it has been argued that close to the Lif- 



4-2. Comparison of different crystalline phases 
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Fig. 4. Comparison of thermodynamic potentials relative to the restored phase for 
different kinds of modulations at T = 0. 

shitz point ID modulations should be favored over higher dimensional ones 
within this type of models [25]. We therefore focus on what happens for 
T = 0, where a GL analysis is unable to give reliable results. 

The results of this comparison are shown in Fig. [H One can clearly see 
that the solitonic solutions ()13p lead to the biggest gain in free energy com- 
pared to all the other cases considered. The two-dimensional "egg-carton" , 
on the other hand, turns out to be energetically disfavored with respect 
to one-dimensional real modulations throughout the whole inhomogeneous 
window, while still being favored over the chiral density wave ansatz ()12p . 

While a thorough discussion will require the study of several kinds of 
2D modulations, these preliminary results seem to indicate that the gain in 
condensation energy is not enough to compensate the greater kinetic energy 
cost due to the presence of a condensate varying in more than one spatial 
dimensions, and thus ID modulations remain the most favored. 



In order to mimic features of confinement, in particular to suppress the 
contribution of free constituent quarks in the confined phase and to in- 
clude gluonic contributions to the pressure, the NJL model can be coupled 
to an effective description of the Polyakov loop [271 l28| I29j . To this 
end, the quarks are minimally coupled to a background gauge field. Fur- 
thermore, a local potential which is essentially constructed to reproduce 
ab-initio results of pure Yang-Mills theory at finite temperature, is added to 
the thermodynamic potential. The resulting model is known as the PNJL 
model, and within the mean-field treatment the traced expectation value of 



5. PNJL and large Nc 
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Fig. 5. PNJL phase diagram varying the number of colors: Nc — 3 (left), Nc = 10 
(center), Nc = 50 (right). The shaded region indicates the inhomogeneous phase. 



the Polyakov loop, £, becomes a new quantity which has to be determined 
self-consistently. 

This extension is also possible when dealing with inhomogeneous phases 
and has been investigated in Ref. pHj under the simplifying assumption 
that the Polyakov loop expectation value i is spatially uniform. It was 
found that the effects of coupling to the Polyakov loop basically amount 
to stretching the phase diagram towards higher temperatures, since the 
effective confinement suppresses the thermal excitation of single quarks. 

On the other hand, the PNJL model has also been used to investigate 
the influence of the number of colors on the phase diagram |30] . Here the 
coupling to the Polyakov loop is crucial, since the NJL mean-field results 
without Polyakov loop are A^c-independent. The study of Ref. [30j was 
restricted to homogeneous phases, and its main focus was to see how the re- 
gion of confined but chirally restored matter develops. Some time ago, this 
was thought to be a possible manifestation of quarkyonic matter, whereas 
according to the present picture chiral symmetry is inhomogeneously bro- 
ken in the quarkyonic phase [8l[9l[T0]. Therefore, since the latter has been 
derived in the large- A'c limit, it is particularly interesting to combine the 
approaches of Refs. [14J and [3D], and to study the behavior of the inhomo- 
geneous phase in the PNJL model at large Nc- 

As the exact implementation of the Polyakov-loop expectation value is 
not unique at arbitrary Nc, we follow Ref. |30j and change the function fmed, 
Eq. (|16p. in the thermodynamic potential into 



fmed = 0(i?,-^)^T(e-(^^-'^)/^ + e-(^^+'^)/^^ 

+ e(/x - Ep) \{fi -Ep) + £T (e-^^"-^"^/^ + e-^^^^"^/^) , (20) 



which basically amounts to a leading-order expansion for small i and is as 
such expected to be accurate at large Nc in the confined phase. For the 
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Fig. 6. Left: Phase diagram with inhomogeneous island and continent. Right: 
Amphtude and wave number for the CDW ansatz at T = using Pauh-ViUars 
regularization. 



shape of the inhomogeneous modulations we consider the one-dimensional 
solitonic ansatz, Eq. p3]) . 

Our results for Nc = 3, 10 and 50 are shown in Fig. [SI By increasing 
the number of colors, the inhomogeneous phase is enlarged and stretches 
towards higher temperatures, approaching the upper limit given by the pure 
glue transition temperature (Tc = 270 MeV in our parametrization of the 
Polyakov loop potential). The transition lines between homogeneous and 
inhomogeneous phases become more and more vertical, assuming a shape 
resembling the expected form for the phase diagram in the large Nc limit 
[31j . The size of the inhomogeneous phase is not dramatically enhanced in 
the ^ direction. In fact, in the PNJL model the Polyakov loop decouples 
from the NJL sector at T = 0. Therefore, since the latter is A'^c-independent 
in mean-field approximation, the transitions at T = are unchanged. 



6. The inhomogeneous continent 

An unexpected feature emerges, when we extend our studies of the phase 
diagram to higher values of As shown in Fig. [6l above some critical value 
of the chemical potential a second inhomogeneous phase appears, which 
seems to extend to arbitrarily high /x with steadily growing amplitude and 
wave number of the spatial modulation. This inhomogeneous "continent" 
is present for all kinds of spatial modulations we considered, and its onset 
is again given by a second-order phase transition. In this case, GL analyses 
reveal that the phase boundary is the same for all inhomogeneous solu- 
tions [23j. Within some parametrizations, especially when one chooses a 
stronger coupling to enforce a larger value of the constituent quark mass in 
vacuum, the continent turns out to be even directly connected to the inho- 
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mogeneous island discussed until now. In fact, indications for this behavior 
are already visible in Ref. |13] . 

Since the continent appears at relatively high chemical potentials and 
with large wave numbers, the first natural interpretation for it is that it is 
an artifact of the regularization. While this might indeed be the case, we 
will argue that the effect is at least non-trivial. In order to achieve some 
better understanding on the observed behavior, it might be worth analyzing 
in detail the mechanisms that within the model lead to the formation of 
an inhomogeneous condensate. Since, as argued above, we expect these 
considerations to be roughly independent of the particular ansatz chosen 
for the spatial modulation, we focus on the simplest possible case, namely 
the CDW (USD. 

In this case, we obtain for the thermodynamic potential 

n-nrest = -NfNc / dE[p{E,M,Q)-prest{E)][f,ac{E,A) + fmed{E)] + — , 

(21) 

where the integral corresponds to the kinetic term i^kin, Eq- (|14p . and the 
M^-term to the condensate term ^cond, Eq. 1^. We have subtracted the free 
energy of the restored phase and explicitly indicated the dependence of the 
different terms on the Pauli-Villars cutoff A, the amplitude M and the wave 
number Q of the CDW. In particular we note that only the function f^ac 
depends on the regularization whereas only the density of states depends 
on the wave number. 

Since the condensate term always disfavors (homogeneous or inhomo- 
geneous) chiral symmetry breaking, a necessary condition for a non-trivial 
phase is that the first term is negative, i.e. the integral must be positive. 
The functions which enter the integrand are displayed in Fig. [7] for several 
cases. As shown in the left panel, the Pauli-Villars regularized vacuum func- 
tion fyac (red solid line) is negative for all energies. Hence, in vacuum, the 
density of states minus its value for the chirally restored solution should be 
negative as well to obtain a positive integrand. Indeed, as seen on the right, 
for Q = the difference p — prest is always negative and its absolute value 
increases with increasing M. Thus, a larger constituent quark mass leads 
to an increase of the free-energy gain in the kinetic term, which is stabilized 
by the condensate term at some optimum value of M. This is the usual 
mechanism for spontaneous chiral symmetry breaking in vacuum. 

In contrast to fvac, the function fmed is positive. Hence, when medium 
effects are included, the sum f^ac + fmed becomes less and less negative 
at small energies until it eventually changes sign. As an example the case 
T = 0, p = 600 MeV is shown in the left panel of Fig. [7] (blue dashed 
line). Since at T = the medium only contributes to energies E < p, 
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Fig. 7. Functions entering the integrand in the thermodynamic potential, Eq. (|2ip . 
Left: The red sohd hne represents the Pauh-Villars regularized function /„ac, cor- 
responding to the Dirac sea contribution, while the blue dashed line indicates the 
sum fyac + fmed at T = and ^ = 600 MeV. The two lines coincide above E = ji, 
since the medium does not contribute beyond that point. Right: Density of states 
p{E,M,Q) relative to the restored case, prest{E) = p{E, 0,0), for homogeneous 
condensates with M = 300 MeV (red dashed hne), M = 100 MeV (green dash- 
dotted line), and for a CDW with M = 100 MeV, Q = 1000 MeV (blue sohd 
line). 



the sum f^ac + fmed remains negative at large energies, but altogether the 
formation of homogeneous condensates, related to a negative function p — 
Prest, becomes progressively disfavored by the medium contributions, thus 
leading to chiral restoration if we restrict ourselves to Q = 0. 

The situation changes, however, when we allow for a CDW with Q > 
2M. In this case p — prest is positive at small energies and changes sign at 
E K Q/2- M2/(2Q) (see Fig. El blue solid Hne on the right). Hence, by 
properly choosing Q, it can be achieved that the factors f^ac + fmed and 
P — Prest have the same sign for all energies, thus more or less optimizing 
the free energy gain in the kinetic part of the thermodynamic potential. As 
fvac + fmed changes sign slightly below E = p at large p, this estimate yields 
Q 7^ 2p, for the favored value. 

These arguments support the statement that the formation of an inho- 
mogeneous chiral condensate is a medium- induced effect [H [TTl [12] . In our 
case it is technically related to the fact that both, p—prest with Q > 2M and 
/medi are positive at small energies. Therefore, since none of these functions 
but only the vacuum term is affected by the regularization, the appearance 
of the continent does not immediately look like a regularization artifact. 

Of course, the dynamical formation of the inhomogeneous condensate 
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eventually depends on the interplay of fuac and fmed, and the explicit ex- 
pression for the Dirac sea term will naturally influence the quantitative 
details. For example, by allowing for a larger cutoff A, the vacuum part 
becomes larger and the second continent is shifted to higher chemical po- 
tentials. However, the same is true for the standard chiral phase transition 
between homogeneous phases, where this is usually not considered to be a 
major problem. 

The existence of the inhomogeneous continent is not restricted to Pauli- 
Villars regularization, but the very same behavior is also present when the 
vacuum term is regularized following the Schwinger proper-time prescrip- 
tion. On the other hand, when the thermodynamic potential is regularized 
by introducing a three- momentum cutoff, it is clear from Eq. (jlip , that a re- 
striction of the in- and outgoing momenta also restricts the size of the wave 
number to which they can couple. Hence, large values of Q are strongly 
disfavored and the second inhomogeneous phase, if it exists at all, cannot 
extend to arbitrary large values of fj,. However, unlike the previous exam- 
ples, this is obviously a regularization effect, since the suppression of large 
Q is directly caused by the cutoff. In fact, for this reason the use of a 
three-momentum cutoff was abandoned in Ref. ^16j in the context of inho- 
mogeneous color super conductor 80 

An example where large values of Q are suppressed in a seemingly more 
physical way is the quark-meson (QM) model, defined by the Lagrangian 



J^QM = ^P{i7^d^-g{a+i-f''Ta7T''))i;-^{d^ad^'a+d^TT^d^'7ra)+U{a,7r) . (22) 

Here a and tt are elementary meson fields and U is a Mexican-hat type meson 
potential, leading to spontaneous chiral symmetry breaking. While NJL and 
QM model are formally quite similar, in the latter the vacuum contribution 
is usually omitted when performing a mean-field treatment. On the other 
hand, for CDW ansatz, which for the QM model reads a{x) = f-j^ cos(Q • x), 
7r^(x) = /7rSin((5 • x), the kinetic term of the mesons yields a contribution 

r\mesons ^ f'^r^'^ /'OQ\ 

"/cm — 2^^^^ ' ^ ' 

which suppresses large values of Q (see Ref. [5j for more details). As a 
result, the inhomogeneous continent is not present in the QM model, as far 
as we can tell numerically. 

* For the same reason the authors of Ref. [TO] introduce a form factor, which restricts 
the quark momenta to the vicinity of the Fermi surface, rather than to smaU absolute 
values. 
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Prom this one might naively conclude that the emergence of the conti- 
nent in the NJL model is due to absence of derivative terms in the La- 
grangian. However, it is well known that in the NJL model the con- 
tribution (j23p is already contained in the vacuum part of the thermody- 
namic potential, which carries a dependence on Q through the energy spec- 
trum O [121 [32]. Indeed, if one expands 

^vac{M, Q) = n^ac{M, 0) + /3„ac,2 + ^vac,4 + ■ ■ ■ (24) 

one can show in a regularization independent way that the so-called spin 
stiffness is given by 

(3vac,2 = \fl , (25) 

giving rise to a term like Eq. (j23p in the thermodynamic potential. Hence, if 
we were allowed to neglect all other terms in Eq. ()24p , we could conclude that 
large values of Q and, thus, the inhomogeneous continent are suppressed in 
the same way as in the QM model. However, in the continent region, Q is 
certainly not small (compared with M or /^). In any case, it is clear that 
the stability against large values of Q cannot be discussed in terms of a 
Taylor expansion for small Q. 

We therefore look at the higher orders in the series in order to get an idea 
how these additional contributions (which become more and more relevant 
at higher values of Q) influence the previous considerations. In particular 
the coefficient of the Q^-term is dimensionless and therefore expected to 
stay finite, even if the cutoff is sent to infinity. For instance, if we employ 
proper-time regularization (see e.g. Ref. |12)). we obtain 

P^acu - -2NfN,e -j^^^ - -^ + 0{^). (26) 

The most important observation is that this term is negative, i.e. it weakens 
the effect of the Q^-term. Although the value of the coefficient is relatively 
small, the Q^-term dominates the Q^-term when Q is of the order of 10 
times /tt. 

In Fig. [8] we compare the full proper-time regularized vacuum thermo- 
dynamic potential with the results of a Taylor expansion to the orders 
and Q"^. In agreement with our considerations above, the Q^-result gets 
reduced by the higher orders. Although the effect is overestimated by the 
Q^-term, this leads to a much more moderate increase of ilvac in the con- 
tinent region. As a consequence, the vacuum contributions do not disfavor 
the formation of a CDW strongly enough in this region to compete with the 
favoring effects of the medium contributions. 

A very systematic investigation of the large-Q behavior of the vacuum 
effective potential has been performed more than 20 years ago in Ref. 
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Fig. 8. Proper-time regularized vacuum thermodynamic potential for a CDW with 
fixed amplitude M — 100 MeV as a function of the wave number Q, using param- 
eters from [12]. The value for Q = has been subtracted. The full result (black 
solid line) is compared with the Taylor series truncated at order (red dashed 
line) and order (blue dotted line). 

Although it turned out that the exact behavior beyond the universal Q^- 
order is strongly regularization dependent, in all cases considered, ^Ivac even- 
tually becomes negative, meaning that even the vacuum is unstable against 
the formation of a CDW with very large Q. While this is clearly unphysical 
and should therefore be ignored, there is no reason why corrections to the 
Q^-terms should not be present at all. Unfortunately, these corrections are 
very model dependent and theoretical input from outside is needed to pin 
them down. In fact, inhomogeneous phases at arbitrary large /i have been 
predicted for the 1 -|- 1-dimensional Gross-Neveu model |19] as well as for 
QCD in the large-A^c limit [6j!|. Therefore, the continent may be not as 
exotic as it appears. 

7. Conclusions 

We presented some of our recent results related to inhomogeneous chiral 
symmetry breaking phases in the NJL model. Thereby one focus was a com- 
parison of different spatial modulations of the chiral condensate, including 
an "egg carton" -like two dimensional ansatz and several one-dimensional 
functions. For all of them we observed the emergence of an inhomogeneous 
"island" around the region where the usual first-order chiral phase transition 
would occur for homogeneous phases. 

A comparison of the free energies at T = seems to support the idea 
that one-dimensional modulations are favored over higher dimensional ones. 



^ For Nc — 3, they are, however, disfavored against color superconductivity [7]. 
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In particular, solitonic solutions [13j inspired by analytical studies of 1+1- 
dimensional models seem to constitute the favored state. While only one 
kind of two-dimensional ansatz has been considered so far, it may very well 
be that this feature is valid in general for all kinds of 2D modulations. In 
fact, it has already been argued in a Ginzburg-Landau study that this is 
the case in proximity of the critical point |25j . 

The inhomogeneous island also exists in the Polyakov-loop extended 
NJL model. In this framework we have studied the effects of varying the 
number of colors on the phase diagram. We found that with increasing Nc 
the inhomogeneous phase is stretched towards higher temperatures so that 
the boundaries between homogeneous and inhomogeneous phases become 
more and more vertical. On the other hand, the size of the inhomogeneous 
phase in the direction is not enhanced dramatically. 

Finally, we investigated the origin of the inhomogeneous "continent" , i.e. 
a second inhomogeneous phase which appears in our calculations at higher 
chemical potential and does not seem to end. We discussed the reliability of 
our model in the region where the continent appears, focusing on effects of 
the regularization process. It turned out that the continent is not the result 
of a trivial cutoff artifact but, unfortunately, regularization dependencies, 
including a known vacuum instability [33) . preclude the model from giving 
definite answers about the phase structure in this regime. 

We thank the organizers for a very interesting workshop. Travel sup- 
port by HIC for FAIR (S.C.) and by the DFG under contract BU 2406/1-1 
(M.B.) is gratefully acknowledged. This work was partially supported by the 
Helmholtz Alliance EMMI, the Helmholtz International Center for FAIR, 
and by the Helmholtz Research School for Quark Matter Studies H-QM. 
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